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The beam-beam interaction is one of the dominant sources of emittance growth and luminos- 
ity Ufetime deterioration. A current carrying wire has been proposed to compensate long-range 
beam-beam effects in the LHC and strong localized long-range beam-beam effects are experimen- 
tally investigated in the RHIC collider. Tune shift, beam transfer function, and beam loss rate are 
measured in dedicated experiments. In this paper, we report on simulations to study the effect of 
beam-wire interactions based on diffusive apertures, beam loss rates, and beam transfer function us- 
ing a parallelized weak-strong beam simulation code (bbsimc). The simulation results are compared 
with measurements performed in RHIC during 2007 and 2008. 



PACS numbers: 29.20.db, 29.27.Bd 



I. INTRODUCTION 

Long-range beam-beam interactions are known to 
cause emittance growth or beam loss in the Tevatron 
[J, Q and are expected to deteriorate beam quality in the 
LHC. A possible remedy to reduce their effects is to in- 
crease the crossing angle. However, it has a side effect of 
reducing the luminosity due to the reduction in geometric 
overlap. Compensation of long-range beam-beam inter- 
actions by applying external electromagnetic forces has 
been proposed for the LHC [3]. At large beam-beam sep- 
aration, the electromagnetic force which a beam exerts 
on individual particles of the other beam is proportional 
to ^, which can be generated and canceled out by the 
magnetic field of a current-carrying wire. Direct-current 
wires were installed in SPS [4], DA^NE [5,], and later 
in RHIC [7] for tests. Results of the SPS wire excitation 
experiments have been reported earlier [6,8^. During the 
KLOE experiment in DA<I>NE, an improvement of beam 
lifetime without luminosity loss due to the compensation 
of parasitic collisions has been observed [5]. Experiments 
in RHIC supplement the SPS and DA^NE tests. The 
beam lifetime in RHIC is typical for a collider storage 
ring and better than in the other two machines. Unlike 
the SPS, two beams circulate and interact in RHIC and 
compensation of a single long-range interaction is pos- 
sible. The energies and other beam conditions are also 
different from the SPS and DA<I>NE. So the results re- 
ported here will result in simulations being benchmarked 
under a wider range of operating conditions. 

Two current carrying wires, one for each beam, have 
been installed between the magnets QS and QA of IPG in 
the RHIC tunnel. Their impact on a beam was measured 
during the physics run with deuteron and gold beams. No 
attempt was made to compensate the beam-beam inter- 
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action since so far only ion beams were available which 
have a beam-beam parameter about a factor of 3 smaller 
than proton beams. However, the experimental results 
help to understand the beam-beam effects because the 
wire force is similar to the long-range beam-beam force 
at large separations. In this report we discuss the results 
of numerical simulations of a wire acting on a beam in 
RHIC using a multi-particle tracking code and compare 
with measurements in RHIC Run 7 and 8 [9]. 

The organization of the paper is as follows: The phys- 
ical models used in the simulation code are described 
in Section [TTl The beam and wire parameters are given 
in Section lllli Section [IV] presents the simulation re- 
sults of the effects of the wire on the beam dynamics and 
the comparison between measurement and simulation re- 
sults. We summarize our results in Section [Vl 



II. MODEL 

In the collider simulation, the two beams moving in 
opposite direction are represented by macroparticles of 
which the charge to mass ratio is that of each beam. The 
number of macroparticles is much less than the bunch in- 
tensity of the beam because even with modern supercom- 
puters it is too time consuming to track 10^^ particles for 
the number of revolutions of interest. The macroparticles 
are generated and loaded with an initial distribution for a 
specific simulation purpose, for example, six-dimensional 
Gaussian distribution for long-term beam evolution. The 
transverse and longitudinal motion of particles is calcu- 
lated by transfer maps which consist of linear and nonlin- 
ear maps. In the simulations, the following nonlinearities 
are included: head-on and long-range beam-beam inter- 
actions, external electromagnetic force by current carry- 
ing wire, multipole errors in the interaction region (IR) 
quadrupole triplets, and sextupoles for chromaticity cor- 
rection. The finite bunch length effect of the beam-beam 
interactions is considered by slicing the beam into several 
sections in the longitudinal direction and by a synchro- 
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beam map fld\. Each shce in a beam interacts with par- 
ticles in the other beam in turn at the collision points. 
In addition, the magnetic field of a finite length of wire 
is calculated using the Biot-Savart law at the particles' 
betatron amplitude. 

The six-dimensional accelerator coordinates x = 

(^x^x\y^y\ z^5^ are applied, where x and y are hor- 
izontal and vertical coordinates, x and y the trajectory 
slopes of each coordinates, z = —cAt the longitudinal 
distance from the synchronous particle, and S = Apz/po 
the momentum deviation from the synchronous energy. 
Coupling between the transverse planes is included in 
the transfer map between elements. We adopt the weak- 
strong model to treat the beam-beam interactions where 
one beam is strong and is not affected by the other beam 
while the other beam is weak and experiences a beam- 
beam force due to the strong beam during the interac- 
tions. The density distribution of the strong beam is 
assumed to be Gaussian. 

It is well known that for a large separation distance 
{':^ a) at parasitic crossings, the strength of long-range 
interactions is inversely proportional to the distance. Its 
effect on a test beam can be compensated by current 
carrying wires which create just the ^ field. The advan- 
tage of such an approach consists of the simplicity of the 
method and the possibility to deal with all multipole or- 
ders at once. For a finite length of a wire embedded in 
the middle of a drift length L and tilted in pitch and yaw 
angles, the transfer map of a wire can be written as 

Dl/2 Te^^Oy, 

where Tg^^Oy represents the tilt of the coordinate system 
by horizontal and vertical angles Ox^Oy to orient the co- 
ordinate system parallel to the wire, is the drift 
map with a length -|, Aik is the wire kick integrated 
over a drift length, and SAx,Ay represents a shift of the 
coordinate axes to make the coordinate systems after and 
before the wire agree. 

At the start of the simulation, the particles in the 
weak beam are distributed over the phase space x = 

^x, x\y, y\ 6^ . The number of particles TV is limited 

by the computational power. In order to make the best 
use of a relatively small number of simulation particles 
compared to the bunch intensity, the initial distribution 
should be optimized. Indeed the initial distribution is 
very important because a proper choice can reduce the 
statistical noise in the physical quantities. The simula- 
tion particles are generated in two steps: 

1. The particle coordinates {x^y^z) of particles can 
be directly generated from the spatial Gaussian 
distribution, p (x, z) = px {x) py (y) pz (z), where 

PC (0 = Pco^^P Since the particle coor- 

dinates are not correlated, one can generate them 



by inverse mapping of each cumulative distribu- 
tion function of horizontal, vertical, and longitudi- 
nal Gaussian distributions, and together with the 
bit-reversed sequence to minimize nonphysical cor- 
relations [12]. 

2. An equilibrium distribution in transverse phase 
space e.g. in the {x^x') plane is p(^x^x^ = 

po exp ^ — 1^ ^— . Since the spatial co- 
ordinate X is determined at the first step, the slope 
X can be obtained from the random variate r of a 
univariable Gaussian, i.e., x = (r — axx) / px- 

Following the above physical model, a beam-beam sim- 
ulation code bbsimc has been developed at FNAL over 
the past few years to study the effects of the machine 
nonlinearities and the beam-beam interactions [13]. If 
required, time dependent effects such as tune modula- 
tion and fiuctuation, beam offset modulation and fiuc- 
tuation, dipole strength fiuctuations to mimic rest-gas 
scattering etc can be included in the model. The code 
is under continuous development with the emphasis be- 
ing on including the important details of an accelerator 
and the ability to reproduce observations in diagnostic 
devices. At present, the code can be used to calculate 
tune footprints, dynamic apertures, beam transfer func- 
tions, frequency diffusion maps, action diffusion coeffi- 
cients, emittance growth and beam lifetime. Calculation 
of the last two quantities over the long time scales of 
interest is time consuming even with modern computer 
technology. In order to run efficiently on a multiprocessor 
system, the resulting model was implemented by using 
parallel libraries which are MPI (interprocessor Message 
Passing Interface standard) [14], state-of-the-art parallel 
solver libraries (Portable, Extensible Toolkit for Scien- 
tific Calculation, PETSc) [H, and HDF5 (Hierarchical 
Data Format) 

III. WIRE AND BEAM PARAMETERS 

The transverse electric field of a round Gaussian bunch 
with transverse rms size a is 

where and are the number of particles and the elec- 
tric charge in the opposite bunch respectively, and r is 
the radial distance from the center of the beam. At a 
small radius r <C cr, the field is proportional to r and 
shifts the betatron tune. This tune shift is characterized 
by the beam-beam parameter ^ = where ro is the 

classical radius of the particle, e = j 0^ is the nor- 
malized emittance. While particles at small r undergo a 
linear tune shift, the particles with r ^ a experience a 
^ force. The long-range effect is nonlinear and may vary 
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Figure 1: (a) Layout of RHIC, and (b) location of the 

long-range beam-beam interaction and wire 
compensators at IPG. Beam-beam collisions occur at 
IP6 (STAR) and IPS (PHOENIX). 



from bunch to bunch if the the bunch pattern is asymmet- 
ric. A current carrying wire generates a magnetic force 
which is (X ^, the same as the long-range beam-beam 
force at large separations. The wire current required to 
compensate a long-range interaction is {I^Ly^) = n^q^c, 
where 1^ is the wire current, and its length. In ad- 
dition to the wire strength, the phase advance from the 
location of the wire to the location of the long-range in- 
teraction should be small for an effective compensation 
of the long-range interaction [17]. In RHIC, a possible 
parasitic collision point is near the separation dipole DX 
magnets near IPG and IPS. There is sufficient space in 
the straight sections between quadrupoles Qs and Q4, for 
the wires. The phase advance between the DX magnet 
and the wire location is 5.7° for (3* = 1 m and is expected 
to provide effective compensation [IS]. 

Figure [T] (a) shows the overall layout of RHIC which 
has two head-on beam-beam collisions at IPG and IPS. 
The Yellow beam revolves in a counter clockwise direc- 
tion, and the Blue beam revolves clockwise. The elements 
and the location of the wire at IPG are shown in Fig. [T] 
(b). Two wires are currently installed in RHIC: one for 
each beam. In these simulations the effects of a current 
carrying wire on a beam have been studied for different 
RHIC runs: gold beam at both injection and store ener- 
gies, and deuteron beam at store energy. The parameters 



of RHIC used in our simulations are shown in Table H 
The length of each wire is 2.5 m. The wire current can be 
varied to a maximum of 50 A. Each wire can be moved 
in the vertical direction over a range of G5 mm. 



IV. EFFECTS OF A WIRE ON A BEAM 



The integrated magnetic fields B from a wire are found 
from [iH 



Bo, 



u — V 



47r 



(3) 



where 1^ is the current of wire compensator, u and v are 

\J (^^)^ + + and \J (%^)^ + + 7/^ respectively, 
and the length of the wire. For the wire kick, the 
the field and the force exerted on the beam are ex 1/r 
where r is the transverse distance from the beam to the 
wire. This force is nonlinear and is expected to have 
a significant impact on the beam quality at sufficiently 
close distances. 

The impact of a wire can be observed in RHIC by mea- 
suring the orbit change, tune shift, the beam transfer 
function and the loss rates. The tune shift is one of the 
fundamental observables and it can be directly verified 
with analytical calculation. However, numerical simu- 
lations allow us to calculate other quantities not easily 
observable but which give valuable insight into the beam 
dynamics and can complement the experiments. These 
numerically calculable quantities include the tune foot- 
print, the frequency diffusion map, the dynamic aperture, 
and the diffusion coefficients to characterize the diffusion 
in action. In this section, we present the simulation re- 
sults of all the above quantities at RHIC injection and 
collision runs, and compare them to measurements where 
possible. 

In the model for dynamics at injection energy we in- 
clude the chromaticity sextupoles, these being the dom- 
inant nonlinearity, and the current carrying wire. It 
should be noted that the persistent current sextupoles 
in the arc dipoles were not included in the model be- 
cause their inclusion would have increased the computa- 
tion time substantially and they have a smaller impact 
on the dynamics compared to the chromaticity correcting 
sextupoles. At collision energy, the beam is strongly in- 
fluenced by both machine nonlinearities and beam-beam 
interactions. In the collision model, we include the chro- 
maticity sextupoles, the nonlinear field errors in the IR 
triplets, and the head-on beam-beam interactions at IPG 
and IPS, together with the wire. There are two variables 
associated with the wire that can be varied: the wire cur- 
rent and the beam-wire separation. In the simulations 
the wire strength is set to the maximum value of 125 Am 
in the RHIC wire compensator so that the effects due to 
the wire may be clearly evident. These results will be 
compared with the data taken at this current during the 
machine experiments of 2007 and 200S. Depending on 
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quantity 


unit 


gold beam 
injection store 


deuteron beam 


energy 


Gev/nucleon 


9.795 


100 


107.396 


bunch intensity (Blue/ Yellow) 


10^ 


0.7/0.7 


1/1 


134/1 


emittance ex,y{9b%) 


mm mrad 


5.8 


18 


17 


at IP6 


m 


(9.93, 10.0) 


(0.97, 0.94) 


(0.87, 0.89) 


{/3x,Py) at wire location 


m 


(119.7, 34.6) 


(1082,392) 


(1194, 393) 


beam-beam parameter ^ 


10-^ 




1.25 


1.46 


revolution frequency 


kHz 


77.8 


78.2 


78.2 


tunes {ux, i^y)'. Blue 




(0.230,0.216) 


(0.220,0.231) 


(0.235,0.225) 


tunes {ux^i^y)'- Yellow 




(0.220, 0.230) 


(0.232, 0.228) 


(0.225,0.235) 


chromaticity {^x^^y) 




(2,2) 


(2,2) 


(2,2) 


max 

vertical separation (dy) 


Am 
m 

mm 
cr 


7-14 


125 
2.5 
3.5 

6-10 


7-12 



Table I: RHIC parameters for gold beam at injection and store, and deuteron beam at store. 



the purpose of the simulation, the beam-wire separation 
is also varied typically in the range of 4-10 a where a is 
the rms beam size at the wire location. The beam-wire 
separation is only in the vertical direction. The Blue wire 
is placed below, the Yellow wire above the beam. The 
simulations reported here will be for the Blue beam. 



A. Tune shift and tune footprint 

A basic check of the simulation is to compare the sim- 
ulated tune shifts with the measured tune shifts and also 
against analytical calculation. The transverse tune shift 
at zero amplitude due to wire kicks is given by 



w 

87r2 (Bp) 0-2 



dl 



(4 



dl] 



(4) 



where dx^y denote the beam- wire distances normalized by 
(Jx,y^ the rms beam sizes at the wire location, IwLw the 
integrated strength of the wire, and Px^y the beta func- 
tions at wire location. At the wire location, the ratio 
of horizontal to vertical beta functions is Px/ Py = 2.8 
for the gold beam at store, and Px/Py = 3.0 for the 
deuteron beam. For IwL^ = 125 Am which is the max- 
imum integrated strength allowed in the RHIC wire in- 
stallation, we get the horizontal tune shifts at 7 <j ver- 
tical separation: Aux = 5.1 x 10~^ for the gold beam, 
and Avx = 5.9 x 10~^ for the deuteron beam at store. 
Without the wire, beam-beam collisions are a major 
source of tune spread and shift. The beam-beam pa- 
rameters for two head-on collisions are ^ = 2.4 x 10~^ 
and ^ = 2.8 X 10~^ for gold and deuteron beams at 
store energy respectively. These tune shifts are calculated 
for zero amplitude particles. For beam- wire separations 
greater than 6 a, the large amplitude particles experi- 
ence a detuning and shift mainly from the wire while the 
tune shift of small amplitude particles is dominated by 
the head-on collisions. 



In our simulations, particle tunes are calculated with a 
Banning filter applied to an FFT of particle coordinates 
found from tracking. Figure[2](a) compares the analytical 
expression in Eq. [Hwith the tune shift from tracking and 
averaged over a Gaussian distribution of particles. The 
changes in the transverse tune as a function of the beam- 
wire vertical separation distance in units of ciy is in good 
agreement with the analytic relation Eq. [H The relative 
difference in vertical tune shift between the simulation 
and theory is ^ 10% dX dy = 7 which stems from our 
usage of the centroid to calculate the tune shift and not 
the zero amplitude particle. The measured tune shifts in 
the recent experiments are compared with the analytical 
values in Fig. EJb). These show that measurements and 
simulations also agree. 

The tune footprint also provides useful information es- 
pecially on the choice of working point and in finding the 
resonances spanned by the beam distribution. Figure [3] 
shows tune footprints from tracking single particles with 
initial amplitudes in the range 0-4 Ux^y for two cases: 
without the wire and with the wire powered at 125 Am 
and 7cr vertical separation from the beam. Comparing 
the two footprints without and with the wire, we observe 
that the horizontal tune shift due to the wire is different 
from the vertical one. The ratio of horizontal to vertical 
tune shifts is, for example, Avxj Avy = 2.8 for the gold 
beam, as shown in Fig. [3] (a), which agrees well with the 
ratio of beta functions at the wire. We observe that the 
width of the footprint is wider at large amplitude due 
to its dependence on the separation distance between a 
target particle and the wire. The wire alone causes the 
horizontal tune spread, for example, at 4 a amplitude to 
be 10 times larger than for the no wire case. Characteris- 
tics of tune footprints of the deuteron beam, seen in Fig. 
[3](b), are similar to that of the gold beam. However, the 
final tunes of the deuteron beam due to the wire is closer 
to the diagonal than the gold beam tunes due to the 
difference of their nominal tunes. Both the 9^^ and 12*^ 
order resonances are spanned by the gold beam while the 
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Figure 2: Plots of tune shift dependence on the wire 
separation distance: (a) simulation and (b) 
measurement [7]. Data set are obtained at gold beam at 
store energy. 



deuteron beam is free from these resonances. The tune 
shifts are slightly larger for the deuteron beam than the 
gold beam because of the difference of beta functions and 
beam sizes. 
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Figure 3: Plots of tune footprints for (a) gold beam and 
(b) deuteron beam at store energy. The no-wire case 
includes sextupoles, IR multipoles, and head-on 
collisions. The wire current is 50 A and wire-beam 
separation distance sets to 7cr. Blue and cyan lines 
stand for 9^^ and 12^^ order resonances respectively. 



B. Frequency diffusion 



We have calculated frequency diffusion maps as an- 
other way to investigate the effects of a current carrying 
wire. The map represents the variation of the betatron 
tunes over two successive sets of the tunes [19]: The vari- 
ation can be quantified by d = log y^Ai^^"^-"Ai^, where 

ISaJx = vi'^ — yi'^ is the horizontal tune variation, in the 
simulations, between the first set of 1024 turns and the 
next set of 1024 turns, and lS.Vy = v^^ — Vy '\ If the 
tunes different from i^y^''^, the par- 

ticle's orbit diffuses. A large tune variation is generally 
an indicator of reduced stability. 

Figure [4] shows the frequency diffusion map of the be- 
tatron tune for the gold beam at collision energy and the 
influence of the wire on the map. The color scale shown 
on the right of each map gives a quantitative measure of 
the diffusion: red color corresponds to larger diffusion, 
and blue color represents less diffusion. In Fig. [H we 
directly see the lines connected to resonances. The wire 
increases the detuning of betatron tune and makes the 
particle motions more chaotic at amplitude beyond 3 a. 



The stability boundary is shrunk further. The diffusion 
map. Fig. [5] (a), of the deuteron beam without the wire 
shows mostly stable motion and only a small region with 
appreciable diffusion. It doesn't show the traces of res- 
onances and looks quite different from that of the gold 
beam. This is to be expected since the resonances span- 
ning the footprint of the deuteron beam are only 12^^ 
order. However, the wire changes the diffusion map sig- 
nificantly, as is evident in Fig. [5] (b). While the region 
of high diffusion for the gold beam is distributed at large 
amplitudes, the regions with large diffusion are observed 
even at small amplitudes with the deuteron beam. The 
red "tongues" even extend to particles in the core, e.g. 
to particles at 0.5 a. This is to be contrasted with the 
tune footprint in Fig [3] which shows that no resonances 
below 12th order are spanned by the beam distribution. 
In later sections we will discuss the correlation of these 
frequency diffusion maps with the dynamic aperture and 
loss rates. 
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Figure 4: Plot of frequency diffusion map of betatron 
tunes (a) without and (b) with wire for gold beam at 
collision energy. Wire strength is 125 Am, and 
wire-beam separation is 7<j. The tune change is 

logarithmically scaled by log J l^v^ + ^^v* 



C. Dynamic Aperture 

Magnet nonlinearities and beam-beam interactions 
limit the dynamic aperture. The dynamic aperture quan- 
tifies detrimental effects of nonlinearities and can be cal- 
culated relatively quickly. The dynamic aperture of an 
accelerator is defined as the largest radial amplitude of 
particles that survive up to a certain time interval; in 
this simulation, we choose 10^ turns for both injection 
and collision which corresponds to 13 seconds circulation 
time in RHIC. It is, for example, only about 10% of the 
RHIC injection period for ions. We examined the depen- 
dence of the aperture on the number of turns by calcu- 
lating it for 10^ — 10^ turns and found that the dynamic 
aperture stays nearly constant after 10^ turns. 

We have seen that a wire has a significant impact on 
the tune shift and tune spread, for example, the tune shift 
due to 125 Am wire strength is 0.005 at 7 a separation. 
The significant increase of the frequency diffusion due 
to the wire implies that the motion of particles becomes 
more chaotic at large amplitudes. Hence we can expect 
that a wire will also have a significant impact on the 
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(b) 

Figure 5: Plot of frequency diffusion map of betatron 
tunes (a) without and (b) with wire for deuteron beam 
at collision energy. Wire strength is 125 Am, and 
wire-beam separation is 7 a. The tune change is 

logarithmically scaled by log x U^v^ + ^^v* 



dynamic aperture. 

At injection energy the simulation model for dynamic 
aperture calculations includes sextupoles, transverse cou- 
pling, as well as the nonlinearity due to the wire. At colli- 
sion energy, the nonlinearities due to the multipoles in the 
IR quadrupoles are added and the head-on beam-beam 
interactions. Particles are distributed uniformly over the 
transverse planes with amplitudes 0-15 a. The linear 
chromaticity is set (^ix^iy) = (+2, +2). Since non-zero 
chromaticity increases the momentum dependent tune 
spread, the off-momentum particles cross more betatron 
resonances during synchrotron oscillations. The plots in 
Fig. [6] show the dynamic aperture for the gold beam at 
injection energy, the gold beam at collision energy and 
the deuteron beam at collision energy respectively. Plots 
(a) and (b) in Fig. [6] show that nonlinearities of the 
IR triplets play a major role in determining the stability 
boundary. When the electromagnetic force due to the 
wire is not present, the boundary at injection energy is 
at 15 cr while this boundary at collision energy falls to 7 a 
primarily due to the IR multipole errors. The effect of the 
head-on beam-beam collisions on the dynamic aperture is 
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negligible over our simulation period since the head-on in- 
teractions do not transport particles to large amplitudes. 
Sextupoles have only a small impact at collision energy 
since their strengths are relatively small compared to the 
IR multipole strengths. The stability boundaries are ap- 
proximately circular, but the wire distorts the boundaries 
near the vertical plane since the beam-wire separation is 
entirely in the vertical plane. The dynamic aperture near 
the vertical plane decreases considerably as the wire ap- 
proaches the beam. The stability boundary along the 
vertical plane is directly proportional to the wire sepa- 
ration while along the horizontal plane, the boundary is 
nearly independent of the separation. The difference in 
the dynamic aperture between the deuteron beam and 
the gold beam at collision energies is approximately 1 a 
when the wire is not present. This can be understood 
from the difference in resonances spanning between both 
beams mentioned in the previous section. With the wire 
powered, the dynamic aperture in the two cases is nearly 
the same. 

We now turn to tune scans of the dynamic aperture. 
One of the key parameters for machine operation is the 
working point. The search for a working point with a 
good beam lifetime and high luminosity is always a ma- 
jor issue. Since the dynamic aperture becomes smaller 
when the betatron tune of beam particles is on or across 
resonances, the tune diagram is often a useful indicator 
of single particle stability. In the simulations, the tune 
scans are performed with increments of = 0.01 in 
the two transverse directions. At each set of the tune 
scan, we load identical distributions in order to avoid the 
uncertainties due to different distribution and beam size 
reduction from the previous tune. 

Figure [3 shows the contour plots of dynamic apertures 
over transverse tunes for two different wire separations 
for gold injection energy. Red indicates low dynamic 
apertures around 3 a while blue indicates high dynamic 
apertures around 11 a. It is found that at all wire sep- 
arations, the largest dynamic apertures are distributed 



along a band parallel to the diagonal, i.e., I'x — h'y 0.02. 
On the other hand, the zone along i^x — — 0-03 has 
the smallest dynamic apertures at all separations. This 
scan indicates that the nominal tune (28.230, 29.216) is 
close to optimal. Furthermore, a sharper drop in dynamic 
aperture is observed near the 5^^ resonances than at other 
resonances as the separation decreases from 8 a to 4 a. 
Figure [8] shows the dynamic apertures over transverse 
tunes for gold collision energy. Here red stands for 5 a 
while blue for 11 a. Since the IR multipoles cause a large 
drop in dynamic aperture at collision energy and the wire 
has a relatively smaller effect, we exclude the IR multi- 
pole errors in order to see the effect of the wire alone. 
As expected, the smallest dynamic aperture is observed 
near 4^^ and 5^^ resonances. At all wire separations, the 
largest dynamic apertures are distributed nearly along 
the diagonal between Ux = 0.21 and Ux = 0.24. The ef- 
fect of the wire on the dynamic aperture of the deuteron 
beam is presented in Fig. [9l It is clear that the re- 
duction of the dynamic aperture is dominant near 4^^ 
resonance. A notable variation is seen near a circular 
band, i.e., -\- Uy 0.21^, when the beam-wire separa- 
tion is small, while the circular band is not detected in 
the gold beam. Indeed, the large stability boundary is 
distributed all over the tune space except the particular 
bands. Since the wire kick at a large separation is equiv- 
alent to long-range beam-beam interaction, these scans 
can be interpreted as showing the effects of a long-range 
interaction at different tunes. 



D. Beam Transfer Functions 

The beam transfer function (BTF) is defined as the 
beam response to a small external longitudinal or trans- 
verse excitation at a given frequency. BTF diagnostics 
are widely employed in modern storage rings due to its 
non-destructive nature. The beam response is observed 
usually in a downstream pickup while either a stripline 
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Figure 7: Tune scan of dynamic aperture at gold 
injection energy: (a) dy = 7.8a and (b) dy = 3.8a. 
Working point (28.230, 29.216) is plotted as a white dot. 

Red indicates low dynamic apertures while blue 
indicates high dynamic apertures. The unit of dynamic 
aperture is a. 
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Figure 8: Tune scan of dynamic aperture at gold 
collision energy: (a) dy = 8 a and (b) dy = 6 a. 
Working point (28.220, 29.231) is plotted as a white dot. 

Red indicates low dynamic apertures while blue 
indicates high dynamic apertures. The unit of dynamic 
aperture is a. 



kicker or RF cavity excites the betatron or synchrotron 
oscillation over the tune spectrum. The fundamental ap- 
plications of BTF are to measure the transverse tune 
and tune distribution by exciting betatron oscillation, 
to analyze the beam stability limits, and to determine 
the impedance characteristics of the chamber wall and 
feedback system [2fl, HH, [22| . Since this is one of the ob- 
servables in RHIC, we calculate the BTF to benchmark 
another output from the code against measurements. 

Figure [10] presents the beam transfer functions 
for deuteron collision energy and nominal tune 
(28.228, 29.225) obtained by applying a sinusoidal driv- 
ing force to a beam in transverse plane and tracking the 
excited particles over 1024 turns at each excitation fre- 
quency of the kicker. The driving frequency is swept from 

(n + 0.21) frev to (n + 0.24) frev in steps of 2 X 10~^/rei;, 

where n is the integer tune and frev the revolution fre- 
quency which are listed in Table [B Through this fre- 
quency scan, we compute the amplitude of beam response 
and its phase. However, as the excitation strength is in- 
creased, it is observed that the response curves are dif- 



ferent according to the direction of the frequency sweep, 
i.e., downward or upward [2^,[23|. The driving amplitude 
should be applied as low as the response can be detected. 
In the simulations, the driving amplitude is chosen as 
10~^ a^' y'^ where a^' y' are the transverse rms slopes. 
Besides, the response profile is affected by the sweeping 
rate. A relaxation time needs therefore to be applied just 
before starting the evaluation of beam transfer function 
at each successive driving frequency, because the direc- 
tional dependence of the response amplitude may be due 
to the remnants of nonlinear oscillations driven by for- 
mer excitations. In order to avoid the uncertainty, we 
reload fresh particles with distribution identical to the 
initial one. 

The response curve reveals two peaks, as shown in Fig. 
[TOl (a) and (b): One peak is close to 0.230 which is the 
horizontal tune, and the other is 0.225 which is the verti- 
cal tune. The amplitude is in arbitrary units and normal- 
ized by the peak amplitude. We observe signs of trans- 
verse coupling in these BTFs, the vertical tune shows 
up as a lower peak in the horizontal BTF and similarly 
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Figure 9: Tune scan of dynamic aperture at deuteron 
collision energy: (a) dy = 8 a and (b) dy = 6 a. 
Working point (28.235, 29.225) is plotted as a white dot. 

Red indicates low dynamic apertures while blue 
indicates high dynamic apertures. The unit of dynamic 
aperture is a. 



the horizontal tune appears in the vertical BTF. It is 
clear that due to resonance with the external force, the 
betatron amplitude of particles grows when the driving 
frequency is close to the betatron frequency. However, 
we observe valleys in the response at Ud = 0.226 and 
Ud = 0.228 for both horizontal and vertical amplitudes. 
They are located at the resonance line of order 9, and 
may stem from the detuning of particles away from the 
resonance. For a sinusoidal driving force, the phase of 
the beam transfer function change its sign through the 
resonance frequency [25]. Since the simulation does not 
account for the phase change due to cable delays between 
the pickup and the network analyzer, the absolute phases 
will not agree between simulations and measurements. 
We can however expect a similar phase variation across 
the resonance. Both simulation and measurement results 
show the flip of the phase angle as shown in Fig. [10] (c) 
and (d). The small differences in the zero crossing phase 
and the differences in peak locations of the amplitudes 
between simulation and measurement are likely due to 
the external impedances and the nonlinearities not in- 



cluded in the simulation. 

The influence of the current carrying wire to the beam 
transfer function is presented through the comparison of 
Fig. m which shows the beam amplitude and phase re- 
sponse obtained without the wire and with the wire pow- 
ered at 50A and at two separations in a deuteron beam 
study. The shift of a peak location of the amplitude in- 
creases as the wire separation decreases while the width 
of the amplitude response widens, as shown in Fig. [TTlfa) 
and as expected from the increase in the size of the tune 
spread. The wire with 10 a and 8 a separations shifts the 
peak location of the horizontal response by 3.4 x 10~^ and 
4.1 X 10~^ respectively, which are equivalent to the tune 
shift of zero amplitude particles. As shown in Fig. [TT] 
(c) , the wire moves the phase proflle leftward commensu- 
rately with the shift of the amplitude. The BTF response 
contains a wealth of information about the beam, a de- 
tailed analysis of measured and simulated BTFs is left to 
a subsequent publication. 



E. Diffusion Coefficients 

We have calculated the beam diffusion due to non- 
linear particle dynamics which includes the nonlinearities 
from the machine itself, the head-on beam-beam interac- 
tions, and the current carrying wire. Growth of parti- 
cle amplitudes may be described by a diffusion in action 
variables. The diffusion coefficients can characterize the 
effects of the nonlinearities present in the accelerator, 
and can be used to find numerical solutions of a diffusion 
equation [26*1. The solutions yield the time evolution of 
the beam density distribution function for a given set of 
machine and beam parameters. This technique enables 
us to follow the beam intensity and emittance growth for 
the duration of a luminosity store, something that is not 
feasible with direct particle tracking. Here we will fo- 
cus only the calculation of the diffusion coefficients and 
compare them with past measurements in RHIC. 

The transverse diffusion coefficients can be calculated 
numerically from 

Ai (auaj) = 1 {{Mai,N) - Ji{auQ)) 

{Jj{aj,N)-Jj{aj,0))), 

where Ji (ai,0) is the initial action at an amplitude a^, 
Ji (ai^N) the action at an amplitude after N turns, 
() the average over simulation particles, and (i, j) are 
the horizontal x or the vertical y coordinates. Equation 
([5]) is averaged over a certain number of turns to elim- 
inate the fiuctuation in action due to the phase space 
structure, e.g. resonance islands. In the simulations, the 
tracking code evaluates the diffusion coefficients in two- 
dimensional action space with the boundary determined 
by the dynamic aperture obtained in the previous section. 
We load the initial particle distribution which is built by 
using 100 particles placed at the same transverse action. 
The tracking is performed for 10^ turns. The diffusion 
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Figure 10: Amplitude of (a) horizontal and (b) vertical beam transfer functions, and phase angle of (c) horizontal 
and (d) vertical BTF for deuteron beam. The amplitude is in arbitrary units, and the phase angle is in degree. The 

wire strength is 125 Am. 



coefficient is averaged every 10^ turns at which the co- 
efficient approaches approximately the asymptotic limit. 
The above process is performed over the transverse action 
space with the boundary typically set at an amplitude of 

7 Cr. 

Figure [12] presents the contour plot of the horizontal 
and vertical components of the diffusion coefficients for 
the gold beam at collision energy. In these plots we ob- 
serve that the diffusion coefficients have angular depen- 
dence. The horizontal diffusion Dxx depends strongly 
upon the horizontal action Jx while its dependence on 
the vertical action is weak. Similarly the vertical diffu- 
sion Dyy depends primarily on the vertical action. The 
diffusion coefficients scale exponentially with initial ac- 
tions, specially at small actions, in both horizontal and 
vertical directions but at very different rates in each di- 
rection. A single exponential fit does not however suffice 
but a combination of at least two different exponentials 
is required to describe the growth of the diffusion coeffi- 
cients from the origin to the dynamic aperture. 

Figure [13] shows the diffusion coefficients for three 
RHIC situations: gold beam at store, deuteron beam 
at store, and proton beam at store. The coefficients are 
plotted at the action J = . / + J2 after averaging them 



at the same action and are compared with the measure- 
ments obtained by fitting the time-dependent loss rate 
after moving a collimator into and out from the beam 
(27| . The loss rate due to the movement of the collima- 
tor can be described by [28(1 

iV(i) (t) = aJl + + ai, 

/ A \ 
iV(2) (t) = aoerfc ( , ^ | + ai, 

where the superscript (1) and (2) stands for the inward 
and outward movement of the collimator, Az the change 
in z due to the collimator movement, R the diffusion coef- 
ficient, and ao and ai are constants. The diffusion coeffi- 
cient R can be obtained directly from the fit of measured 
loss rates. In Fig. [131 the vertical axis is a logarith- 
mic scale. It should be noted that dependence of dif- 
fusion coefficients on the initial action is exponential at 
small amplitudes and power law like at larger amplitudes. 
However, since the measured coefficients are fitted by a 
power law, i.e. D ^ J^, they agree with simulations only 
at large actions. Since the collimators were not moved 
into the beam core, the diffusion coefficients were not 
measured at small actions. Conversely it is difficult to 
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Figure 11: Amplitude of (a) horizontal and (b) vertical beam transfer functions, and phase angle of (c) horizontal 
and (d) vertical BTF for deuteron beam. The amplitude is in arbitrary units, and the phase angle is in degree. The 
wire current is set 50A. The separation is in unit of rms beam size at wire location. 



calculate the coefficients in simulations at large action 
because some of the particles are lost quickly. The ef- 
fects of wire on diffusion coefficients are considerable at 
all action amplitudes, as shown in Fig. [131 For exam- 
ple, the diffusion coefficient at 3 <j amplitude becomes 
20 times larger when the wire with strength 125 Am and 
separation dy = 8 a is applied. First we observe that 
the relative increase of diffusion coefficients at below 3 a 
amplitude for the deuteron beam is higher than that for 
the gold beam. In general the diffusion process depends 
on the particle motion in phase space and the resonance 
structure due to the nonlinearities. The differences in 
tunes change the cross-talk between the different non- 
linearities as well as the resonance driving terms which 
are likely to be responsible for the differences in diffu- 
sion. The enhanced diffusion at near 3 a amplitude for 
the deuteron beam leads to significant increase of parti- 
cle loss under the simulation conditions (see the following 
section). 



F. Beam losses as a function of wire-beam 
separation 

The beam lifetime in RHIC is determined by the in- 
elastic interactions and beam-beam interactions when the 
beams are in collision as well as nonlinearities of machine 
elements, intra-beam scattering and residual gas scatter- 
ing. In this section we will focus on the impact of the 
wire on the beam loss rates as the beam-wire separation 
is changed. In this study, the initial beam distribution is 
a hollow Gaussian in transverse phase space and a normal 
Gaussian in longitudinal phase space. The initial trans- 
verse beam sizes are obtained from the RHIC optics and 
typical initial emittances while the bunch length is taken 
from measured values. 

The tracking is done with 5 x 10^ macroparticles, and 
carried out over 10^ turns for each wire separation. The 
loss rates are estimated from the asymptotic limit by ex- 
trapolating the simulated loss rate from 10^ turns to in- 
finity as shown in Fig. [TH It is observed that in the 
beginning of the simulation, the loss rate decreases expo- 
nentially rapidly and then approaches a constant rate at 
later times. We apply, therefore, an exponential decay fit, 
i.e., ae~^^^ + 6, to the loss rate data between 4 x 10^ and 
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(a) 




(b) 

Figure 12: Contour plot of (a) the horizontal diffusion 
coefficient D^x and (b) the vertical diffusion coefficient 
Dyy for gold beam at collision energy. The coefficients 
are calculated at the Blue ring. The color assignment is 
logarithmically scaled in the plots. 



1 X 10^ turns. The asymptotic loss rate is the parameter 
h. 

Figure [15] plots the asymptotic beam loss rate due to 
the wire as a function of beam-wire separation distance 
for the case of Table H For the gold injection energy 
study, we include the chromaticity correcting sextupoles, 
chromaticity set to two units, and the wire. It can be seen 
in Fig. [6] (a) that the electromagnetic force of the wire 
decreases the dynamic aperture significantly compared to 
the case when the wire is not present. The particle loss 
rate of the beam shows a sharp increase wire separations 
smaller than 8 cr. 

For the collision energy simulations, we include the 
nonlinear field errors in the triplets, the head-on beam- 
beam interactions, and the nonlinearities in the injec- 
tion energy study. The onset of beam losses, seen in 
Fig. [m (b) and (c), is observed at 8 a and 9 a for 
gold and deuteron beams respectively. In all three cases, 
the threshold separation for the onset of sharp losses ob- 
served in the measurements and simulations agree to bet- 
ter than 1 cr. It is also significant that the simulated loss 
rates at 7 and 8 a separation for the gold beam and 8 
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Figure 13: Plot of diffusion coefficients of (a) gold, (b) 
deuteron, and (c) proton stores of RHIC. The 
coefficients are calculated at the Blue ring. The 
wire-beam separation distance is 8 a at a wire location. 
The coefficients were measured and fitted in stores in 
previous years [27]. 



and 9 a for the deuteron beam are very close to the mea- 
sured loss rates. At fixed separation, the wire causes 
a much higher beam loss with the deuteron beam than 
with the gold beam. For example, the loss-rate for the 
gold beam at a 8 a separation is about 10 %/hr while 
for the deuteron beam the loss rate is about an order of 
magnitude higher both in measurements and simulation. 
This difference is not reflected in the dynamic apertures, 
shown in Figure [H which are about the same for the two 
cases at the same beam-wire separation. Nor is this dif- 
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Figure 14: Plot of the variation of loss rate versus the 
number of tracking turns. The loss rate is estimated 
from the asymptotic value of the exponential decay fit, 
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Table II: Dependence of loss rate on the beam-wire 
separation: The loss rates are fitted to r ex d^, where dy 
is the separation, and the power a is listed on the table. 
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ference correlated to the tune footprint and the resonance 
lines, seen in Fig. [3l where the footprint for the deuteron 
beam is free of resonance lines lower than the 12^^ while 
the gold footprint is spanned by the 9^^ and 12^^ order 
resonances. However the frequency diffusion maps with 
the wire, seen in Fig. [4] and Fig. [5l show greater dif- 
fusion in the deuteron case than in the gold case. This 
correlation of the frequency diffusion maps with the loss 
rates observed here deserves to be studied more deeply 
and tested for validity in other accelerators. The action 
diffusion seen in Fig. [13] is also larger in the deuteron 
beam than the gold beam at small amplitudes by one to 
two orders of magnitude. Thus both frequency and ac- 
tion diffusion seem to be better correlated with loss rates 
than the traditional short term indicators like footprints 
and dynamic aperture. 

Changing the beam-wire separation changes several 
parameters including the tunes, the tune spread, the res- 
onance driving terms etc. The wire separation alone can- 
not describe the change in the dynamics that infiuences 
the loss rates. Nonetheless the loss rates seen in Fig. [15] 
can be fitted to a power law in the separation. The re- 
sults for the three cases for the measurements and simu- 
lations are shown in Table HI! The simulations have stud- 
ied smaller separations than would be practical in RHIC - 
the large loss rates seen at the smallest separation would 
have quenched the machine. However, the simulations 
show the steep climb in the loss rate beyond a thresh- 
old separation and consequently will have a higher power 
law behavior. The power law at injection is somewhat 
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Figure 15: Comparison of the simulated beam loss rates 
with the measured as a function of separations, (a) gold 
beam at injection energy, (b) gold beam at collision 
energy, (c) deuteron beam at collision energy. Wire 
strength is 125 Am. 



higher than at collision but the power laws for gold and 
deuteron beams are very close despite the large difference 
in loss rates. We expect that the power laws depends on 
the details of the machine and not to be universal. In 
fact different power laws have been reported for the SPS 
and the Tevatron |i,[2i,[33|. 
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V. SUMMARY 

In order to study the effects of the machine nonhnear- 
ities and the beam-beam interactions, including strong 
locahzed long-range beam-beam interactions, we have de- 
veloped a six-dimensional weak-strong code bbsimc. Ma- 
chine nonlinearities, beam-beam collisions and the field 
of a current carrying wire can be included in the model. 

We have studied the effects of the wire on the beam 
dynamics in three different cases: gold beam at injection, 
gold beam at store, and deuteron beam at store. Results 
show that the betatron tune change due to the wire is 
well tracked by the simulation, and that the stability of 
particle motion is strongly infiuenced by the wire which 
causes a significant increase in tune spread and diffusion 
for both gold and deuteron beam. Tune diffusion with the 
gold beam and the wire is appreciable only at amplitudes 
larger than 2 a while the tune diffusion with the deuteron 
beam and the wire is significantly larger extending to the 
beam core at amplitudes down to 0.5 a. 

At injection energy, the dynamic aperture is largely 
determined by the wire, while the wire and the IR mul- 
tipoles have a major impact on the stability boundary at 
collision. With the beam-wire separation only along the 
vertical axis, it is found that the stability boundary near 
the vertical axis is linearly proportional to the beam- wire 
separation. We observe, from the tune scan of dynamic 
aperture, that at injection the largest and the smallest 
dynamic apertures are distributed along a band parallel 
to the zone along Vx — J^y — 0.02 and Vx — J^y — 0.03 respec- 
tively at all wire separations while at collision the largest 
dynamic apertures are distributed nearly along the diag- 
onal between Vx = 0.21 and Ux = 0.24. By modeling the 
stripline kicker we obtain the amplitude and phase angle 
of the transverse BTF. Results show that the betatron 
tune is well identified by the simulation, and that the 
wire changes characteristics of the beam response signif- 
icantly in the horizontal plane, similar to the result for 
the tune shift. 

In simulations the wire enhances the diffusion. This ef- 



fects is particularly pronounced in the plane of the wire, 
and at amplitudes of about 3 <j, where the wire can in- 
crease the diffusion rate by more than an order of mag- 
nitude. The action diffusion for the deuteron beam is 
larger than for the gold beam, similar to the result for 
frequency diffusion. Simulations of the beam loss rate 
when the wire is present are in good agreement with the 
experimental observations. The threshold separation at 
which there is a sharp rise in the loss rates agree to bet- 
ter than 1 a. The loss rates with the deuteron beam 
are found in measurements and simulations to be nearly 
an order of magnitude higher than with the gold beam. 
Comparisons of tune footprints and dynamic apertures 
of these two cases show no indication that the loss rates 
may be higher with the deuteron beam and in fact the 
tune footprint with the deuteron beam is free of low order 
resonances. However frequency diffusion and action dif- 
fusion with the deuteron beam are substantially higher 
than with the gold beam implying that these diffusion 
measures may be better indicators of loss rates. 

In this paper simulation results of the beam-wire in- 
teractions were compared with measurements taken when 
the long-range beam-beam interactions were not present. 
It is expected that compensation of parasitic interactions 
using the wires will be tested during the RHIC run in 
2009. We plan to compare simulations with these forth- 
coming experimental results and also to determine the 
effectiveness of the compensation of long-range interac- 
tions in the LHC. 
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